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INTRODUCTION 


The  question  of  whether  or  not  electroriagnetic  cyclotron  harmonic 
waves  can  be  excited  in  the  magnetosphere  near  qeostationary  altitude, 
was  raised  upon  examination  of  the  data  from  the  VLF  broadband  measure¬ 
ments  on  the  SCATHA  satellite  during  a  severe  spacecraft  charginq 
event  (Donatelli  et  al.,  1983).  During  this  event  strong  emission 
bands  were  detected  at  or  near  the  first  two  harmonics  of  the  2  kHz 
electron  gyrofrequency.  These  emission  bands  appeared  on  both  the 
electric  and  magnetic  wave  field  detectors.  While  the  SCATHA  electron 
gun  was  in  operation  at  varyinq  current/energy  levels  attempting  to 
discharge  the  satellite,  there  were  frequency  shifts  in  the  emission 
band  and  consistent  maqnetic  field  components.  This  observation  prorp- 
ted  the  present  theoretical  study  to  determine  whether  electromagnetic 
electron  cyclotron  harmonic  waves  could  be  generated  in  a  plasma  repre¬ 
sentative  of  local  conditions. 

Theoretical  and  observational  studies  to  date  have  explained  obser¬ 
vations  analogous  to  those  of  the  SCATHA  data  in  terms  of  electrostatic 
emissions.  Using  the  Harris  dispersion  relation,  Fredricks  (1971) 
showed  that  non-resonant  electrostatic  waves  could  be  excited  in  a 
plasma  in  which  the  electron  perpendicul ar  velocity  distribution  has  a 
narrow  reqion  of  positive  slope,  i.e.  3f0 /8vj_  >  0.  The  restriction 
on  the  velocity  distribution  was  relaxed  by  Younq  et  al .  (1973)  to 
the  requirement  that  a  distribution  of  cold  and  warm  electrons  have  a 
velocity  distribution  that  is  non-monotonic  in  V ±  .  Thev  shewed 
that  this  provides  sufficient  free  energy  for  drivinq  instabilities. 
The  ratio  of  cold  to  hot  plasma  density,  above  which  instability  will 
not  occur,  was  established  by  Ashour-Abdall a  et  al.  (197b).  The  non- 
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convective  and  convective  nature  of  the  instability  is  shewn  to  be 

controlled  by  the  temperature  ratio  of  the  cold  to  hot  electron  popula¬ 
tions  ( Ashour-Abdal la  and  Kennel,  1978).  FI ectrostatic  upper  hybrid 

waves  observed  on  ISF.F-1  have  also  been  explained  by  this  theory  (Kurth, 
et  al.  1979).  Sentnan  et  al.  (1979)  show  examples  of  two-component, 
electron  distribution  functions  constructed  from  ISEE-1  observations 
of  1  ow  energy  maqnetospheric  electrons  that  occur  simultaneous  with 
the  detection  of  electrostatic  emissions  between  the  electron  cyclotron 
frequency  and  the  upper  hybrid  resonance. 

Flectromaqnetic  emissions  within  this  frequency  ranqe,  with  the 
exception  of  the  whistler  mode,  have  not  been  reported  prior  to  the 

SCATHA  observations.  Ohnuma  et  al.  (1981)  have  shown  that  electromag¬ 
netic  cyclotron  harmonic  waves  may  he  generated  in  a  dense  plasma; 
i.e.  one  in  which  wpe  >>  QCP,  where  oipe  and  nce  are  the 
electron  plasma  and  gyro-  frequencies,  respectively.  Although  their 
results  were  applied  to  laboratory  plasmas,  this  condition  may  prevail 
in  the  vicinity  of  SCATHA  during  beam  operations. 

Here  the  problem  will  be  approached  in  a  manner  suggested  by  the 
work  of  Tataronis  and  Crawford  (1970).  They  used  the  quasi -static  ap¬ 
proximation  in  deriving  a  dispersion  relation  for  a  warm  maqnetopl  asna 
and  proceeded  to  examine  under  what  conditions  unstable  electrostatic 
cyclotron  harmonic  waves  ( Re rnstein  nodes)  exist.  They  looked  at  several 
types  of  distribution  functions  and  examined  ranqes  of  plasma  parameters 
for  finding  unstable  modes.  They  show  that  if  the  analytic  solution 
to  the  dispersion  relation  undulates  about  zero  for  a  given  distribution 
function,  non-convective  mode  couplfnq  instabilities  occur.  These  occur 
between  specific  ranges  of  the  two  ratios,  the  plasma  frequency/electron 
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qyrofrequency  and  the  qyrorarlius/wavelenqth,  for  each  harmonic  pair. 
For  two  component  distributions,  the  density  and  velocity  ratios  of 
the  two  components  must  also  be  considered.  Similar  solutions  are 
souqht  here  for  electromaqnetic,  cyclotron  harmonic  waves  (extraordinary 
and  ordinary  modes). 

In  the  next  section  the  electromaqnetic  dispersion  relation  for  a 
warm,  uniformly  maqnetized,  homoqeneous  plasma  is  examined.  Analytical 
solutions  will  he  presented  for  the  three  perpendicul arly  propaqatinq 
modes,  the  Bernstein,  extraordinary,  and  ordinary  modes,  usinq  a  two 
component  distribution  function.  The  ambient  plasma  is  represented  as 
a  Maxwellian;  the  electron  beam  is  represented  as  a  delta-function  in 
velocity  space.  Sample  calculations  show  that  instabilities  may  exist 
for  both  the  Rernstein  and  the  extraordinary  modes.  The  frequencies 
and  wave  numbers  of  the  excited  modes  vary  with  the  beam-to-anbient- 
plasma  density  and  velocity  ratios. 

THEORY 

The  dispersion  relation  is  derived  usinq  the  linearized  Vlasov 
equation  and  Maxwell's  equations  (in  cqs  units): 

[3/3t  +  v  •  7  +  ((v  x  B0)/c)  •  Vy]fji  * 

-  (qi/nj)[Ej  +  (v.  x  ii)/c]  •  Vyfj0  (l) 

_vx  Ej  =  -c"1  aBj/at 

V  x  Bj  =  c-1  aFj/at  +  (Airq j n j ^/c )  /  vdyf  ^ 

V  •  Ej  *  Anq^nji  /  dvf i i 

_V  •  Rj  =  0 
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The  equilibrium  electron  velocity  distribution  is  fjo(vi>vu) 
where  l  and  11  refer  to  components  perpendicular  and  parallel  to  the 
external  magnetic  field,  R0  =  Roz.  There  is  assumed  to  be  no  external 
electric  field,  F0  =  f).  The  perturbation  electric  and  magnetic 
fields  are  Ej  and  Bj ,  respectively.  The  perturbed  particle  distribution 
and  density  are  fji  and  n-j] ,  respectively,  where  particle  species 
is  indicated  by  the  subscript  ,i.  The  charge  and  mass  for  each  particle 
species  are  qj  and  n-j,  respectively;  c  is  the  velocity  of  liqht.  in  a 
vacuum. 

The  plasma  is  assumed  to  be  infinite,  spatially  homogeneous  and 
uniformly  naqnetized,  since  variations  in  time  and  space  nay  be  neglected. 
Equations  (1)  and  (?)  are  solved  by  introducing  a  Fourier  transform  in 
space,  a  Laplace  transform  in  time,  and  integrating  along  unperturbed 
particle  orbits.  This  leads  to  the  nine  element  dielectric  tensor,  from 
Krall  and  Trivelpiece  (1973): 


nxx  nxy  ^xz 

*x 

n  yx  nxy  nyz 

nzx  nzy  r)zz 

where  the  elements  of  the  dielectric  coefficient  are  (lm(u>)>0): 
0XX  =  l  -  (cW/u?)  -  (2u/u))  ^  lN  u>p/  II  xjVnV^i2 
!1xy  =  “  (^lfi/w)  cj).N  10pj2  U  Xj'W 
\z  =  (C^jk/u*)  -  (?»/,o)  ).j).N  u>p/  II  A^J^V  ]|NQj/ki 


Dyx  3  -  11  xy 

nyy  =  1  -[c2(k/-kL2)/J]  -  (2n/ w)  u>p/  II  Xj(.)N,)S? 
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Oyz  =  (2  iri/  a))  (i)pi-2  II  Aj*JN  jn‘  V1VH 

Dzx  =  (c2kjkj/u2)  -  (2ir / a>)  upj^  I  I  Vfl^j^kl 

nzy  =  -  (?"i/u))  wpj-  II  Xi>VVv«vl 

0Z2  =  1  -  (chL2/J)  -  (2u/a))  Ijln  Wp/  II  AjJ^v/ 

where: 

II  =  dvB/0  dV  i2vi/(Nn,i  +  k  II V II  •  w 

Xj  H  [l-(k|1Vfl/a))](3f1-0/3vi2)  +  (kjVn/u)) 

(  3f  jo/3v  11^ ) 

A-j  =  (NOj/uJCOf^/Sv/)  -  (3fj0/3v„2)] 

+  (8fj0/3v|Z) 

For  simplicity  in  the  above  derivation,  the  wave  vector  k_  is 
defined  as  follows: 

i  =  Jit  +  in  *  ki£  *  kii* 

The  definition  of  other  terms  are: 

a)  =  wave  frequency 
'jpj  =  pi  a sma  frequency  =  (A^q 
fj-j  =  particle  qyrofrequency  =  q-jPg/njc 
»)N  =  Nth  order  Ressel  function  of  the  first  kind 
kivi/n.i=  arqument  of  the  Bessel  function 
Jpi ’  =  derivative  of  Juj  with  respect  to  the  arqument. 

There  is  no  separation  of  electrostatic  and  electronaqnetic  modes 
in  the  dispersion  tensor  (3).  However,  by  considering  parallel  propa¬ 
gation  (kL=  0)  and  perpendicular  propagation  (k8  =  0)  separately, 
it  is  possible  to  simplify  the  dispersion  tensor  and  examine  some 
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characteristic  wave  nodes.  Setting  kt  =  n  in  (3)  singles  out  the 
nodes  that  propagate  parallel  to  the  equilibrium  magnetic  field. 
These  nodes  are: 

(1)  The  longitudinal,  electrostatic,  Landau  damped,  Lanqimir  waves 
and  ion  acoustic  waves. 

(?)  The  weakl.v  damped  transverse  electromagnetic  Alfven  waves  at 
frequencies  below  the  ion  cyclotron  frequency;  the  Whistler  waves  at 
y  ~  (1/2 )ft  and  the  cyclotron  waves  which  occur  at  u  =  ft.  The  cyclotron 
waves  are  stronqly  damped  for  w  =kt|C,  which  is  also  the  hand  in  which 
spontaneous  emission  of  radiation  occurs. 

Ry  taking  k  „  =  0,  the  characteri st i c  waves  that  propagate  perpen¬ 
dicular  to  the  equilibrium  maqnetic  field  are  singled  out.  These 
modes  are: 

(1)  the  almost  purely  longitudinal,  nearly  electrostatic 
Rernstein  modes  that  correspond  to  the  electrostatic  waves  discussed  at 
the  beginning  of  this  section. 

(?)  The  transverse,  electromagnetic,  cyclotron-harmonic  waves 
which  include  the  extraordinary  mode  with  [  i  B0  and  the  ordinary 
mode  with  £  ii  R_0*  These  modes  are  normal  modes  of  a  high  density, 
maqneti zed  plasma  (^  >>  fte?),  (Ohnuma  et  al.,  1981)  and  are 
excited  near  ion  and  electron  cyclotron  harmonics.  The  SCATHA  observa¬ 
tions  are  within  the  range  of  these  emissions.  It  will  be  determined 
if  these  modes  are  unstable  within  the  local  beam-backg round  plasma. 

For  k H  =  0,  the  dispersion  tensor  (3)  reduces  to: 


IJxx  ^x.y  0 

-D  Xy  Oyy  f) 

0  0  Dz  z 


n 


(4) 
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At  the  frequencies  of  oscillation  considered  here,  ions  nay 
be  regarded  as  providing  a  charqe  neutral  background.  The  approximation 
<<  kzcz  is  also  valid  since: 

<^2/kzc2  «  u)p2vh2/!D2c2  <  idp^v^/aV  s  0.1 


where:  vp  =  beam  velocity. 

Since  the  0xy  terns  are  second  order  in  (OpZ/k2cZ,  they  may  be  neqlected. 
Then  (4)  may  be  approximately  solved  for  the  three  eigenmodes: 


(1) 

nxx  =  f,: 

Bernstein  node 

(?) 

0yv  =  0: 

extraordinary  mode 

(3) 

nzz  =  0: 

ordinary  mode 

where  the 

elements 

of  the  dispersion  matrix 

reduce  to: 

nxx 

=  1  - 

(4n<dp2/k2 

)  l  N2nz/w(No  -  u)  J  dv, 

{  dv 

1  viV  »Vavl? 

(5) 

nvy 

=  1  - 

(kZcV) 

■  4irt,)b?  1  (w(Nn  -  w))'1 

/  dvn 

/  dv 

1  v,3^' 

)?-  af0/3viz 

(«) 

nzz 

«  1  - 

(k2cZ/wZ) 

-  4l,lj)p?  \  (u)(Na  -  u)))'1 

I  dvn 

/  dv 

l  vl  '? 

JN2[Nn/(d  (3f0/3v/ 

-  3f 

o'*'?)  - 

»f0/»V,Z)3 

(7) 

Onstable  solutions  are  now  sought  for  each  of  the  eigenmodes 
usinq  the  SCATHA  data  to  evaluate  the  necessary  parameters. 

ANALYTICAL  SOLUTIONS 

The  SCATHA  observations  presented  in  Donatelli  et  al.  (1083)  for  ?A 
April  1 ')  7  ,  indicate  that  beatn-  injected  electrons  may  create  a  dense 


pi  a sna  where  electromagnetic  cyclotron  harmonic  instabilities  are  gener¬ 
ated.  At  the  tine  of  interest  the  maqnetnspheric  plasma  was  made  up  of 
two  populations:  a  low-energy  component  with  a  temperature  about  3 on  eV, 
and  a  high-energy  component  about  2F>  keV  (Mullen  et  al.,  1981).  The 
ejection  energy  of  the  beam  from  the  electron  qun  was  either  GO  eV  or 
ISO  eV.  After  passing  through  the  satellite  sheath  the  bear,  electrons  had 
energies  of  1-3  keV  due  to  acceleration  throuqh  the  vehicle  potential. 
Emission  bands  were  detected  consistently  at  or  near  the  first  and 
second  harmonic  of  the  electron  gy rof requency ,  £2,  which  was  about  2  kHz. 
From  the  data,  bean-to-backq round  velocity  ratios  are  estimated.  A  ratio 
of  2-3  is  reasonable  for  the  "artificial"  beam  (electrons  from  the  bear 
systems  on  SCATHA),  and  a  ratio  of  9-10  for  the  natural  beam  (injec¬ 
tions  of  high  energy  naqnetosperic  electrons).  The  density  ratios  are 
assumed  to  be  greater  than  one,  both  for  the  natural  beam  (Mullen  et 
al.,  1981)  and  the  "artifical"  beam.  Although  the  current  and  radius 
or"  the  artificial  beam  are  known,  the  effective  ueam  density  is  not, 
since  the  el ectrostatic  forces  between  bean  electrons  contribute  to 
rapid  spreading  as  does  the  external  magnetic  field  (Gendrin,  1973). 
Furthermore,  at  these  low  emission  energies  the  electron  bean  cannot  be 
highly  focussed. 

The  dimensionless  variables  to  be  used  in  these  solutions  are  de¬ 
fined  as  fol 1 ows: 

s?  =  2KTk?/in«2  =  vTk?/£2? 
xb  5  vib/yT 

where: 

v tb  =  vb  sin  <j>b; 

^  =  pitch  angle  of  the  electron  bean,  15°  <  <j>  «  165°. 


and: 


sxfc  =  kvxb/Q  =  (kvb/fl)sin  <j,b 


0.26  (kvb/£l)  <  sxb  <  kvb/a 

where  sxb  is  the  arqunent  of  the  Ressel  function  in  equations  (5),  (fi)  , 
and  (7). 

The  ambient  density  is  between  f).  5  and  1.0  cnr^.  If  the  density  is 
set  at  1.0  err 3,  the  value  of  the  followinq  non-dimensional ized 
parameters  are: 

(up/ii)2  -  20 

(o^/kc)2  -  0.1 

The  distribution  function  is  approximated  as  follows: 

f0  =  np(n/?ffKT)3/2exp[-(n/2KT)  (v  2  +  v^)]  + 

(nb/2Ttyib)  «(v,  -  v,b)  <5(vl  -  v ib )  (8) 


where: 

T  =  temperature  (°K) 

K  =  Roltznann  constant 

nD  =  1  cw  energy  ambient  electron  density 

nb  =  density  of  election  heam 

This  distribution  function  is  a  Maxwellian  combined  with  a  rinq 
distribution  in  velocity  space  and  includinq  motion  parallel  to  the 
magnetic  field.  The  Maxwellian  represents  the  ambient  plasma,  and  the 
delta-function  represents  the  beam,  with  the  ring  distribution  descri¬ 
bing  the  portion  of  the  monoenerqeti c  bean  electrons  movinq  in  the 
plane  perpendicular  to  the  magnetic  field,  uniformly  distributed  in 
gyrophase  angle.  Tataronis  and  Crawford  (1970)  conducted  a  numerical 
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study  of  t.ht-'  propaqat  i on  charuc t  or i st. ics  of  perpend  it  til  arly  propagating 
electrostatic  waves  in  a  combined  rinq-Maxwel 1 ian  di  strihutior,.  They 
state  that  instability  occurs  if  the  analytic  function  for  the  disper¬ 
sion  relation  is  undulatory  about  2.ero.  Their  results  show  that  node 
couplinq  is  a  feature  of  the  ring  distribution  leadinq  to  stronq  non- 
convective  instability.  Combining  the  rinq  and  Maxwellian  distribu¬ 
tions  leads  to  non-convective  instability  with  higher  qrowth  rate  and 
lower  instability  threshhold  than  for  the  ring  distribution  alone. 
The  electromagnetic  dispersion  relation  nay  be  solved  for  each  of  the 
three  eigenmodes  of  perpendi cul arly  propagatinq  waves.  Solutions  to 
the  dispersion  relation  nay  be  found  as  functions  of  xp  and  np/np, 
usinq  precedinq  definitions. 

A.  The  Bernstein  Mode 

The  dispersion  relation  for  the  Bernstein  node,  obtained  by  substi¬ 
tuting  equation  (8)  into  equation  (5)  and  integrating,  is  the  following: 

Hxx  =  1  -  (4.^/il2)  \  N2a2ANx/(u>2  -  N2fi2)  =0  (9) 

where: 

aNx  =  s~2exp(-s2/2)IN(s2/2)  +  (nb/npsxb)JNJN' 

lf,  is  the  modified  Bessel  function  of  the  first  kind  with  arqument  s2/2. 
The  prime  denotes  the  derivative  of  the  Bessel  function  with  respect  to 
the  argument.  \  is  now  the  summation  from  N=1  t.o  «.  Since  Aryj x 
is  undulatory  about  zero,  unstable  solutions  are  anticipated  (Tataronis 
and  Crawford,  1970).  In  Appendix  A  it  is  shown  that  in  usinq  a  three- 
term  approximation  to  equation  (9),  unstable  solutions  are  found  for: 

nb/np  =  ?  and  xp  -  2 
such  that: 
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iV  fi  =  1.7  +  0.  Pi 


(l(') 


and: 

nb/np  =  4.5  and  x^,  =  in 

such  that: 

•Vfi  *  1.8  +  l.li  (11) 

The  Bernstein  Mode  is  nearly  a  pure  electrostatic  mode.  The 
solutions  are  equivalent  to  those  obtained  by  Tataronis  and  Crawford 
(197(1)  usinq  the  electrostatic  approximation  to  the  dispersion  relation. 
They  found  unstable  solutions  in  this  frequency  ranqe  with  the  qrov/th 
rate,  y,  a  finite  fraction  of  the  real  part  of  the  frequency,  «r.  In 
equation  (10),  with  u>r  =  1.7  fl,  y  =  0.5u>r;  for  equation  (11), 
y  =  0.6^.  In  both  equations  (9)  and  (10)  wr  is  in  the  range 
of  emissions  detected  by  the  SCATHA  broadband  receiver  (Donate!  li  et. 
al.,  1 983) .  However,  this  node  is  not  expected  to  have  the  observed 
naqnetic  field  component,  the  presence  of  a  maqnetic  component  requires 
the  existence  of  electromaqnetic  extraordinary  and/or  ordinary  modes. 

D.  The  Extraordinary  Mode 

The  dispersion  relation  for  the  extraordinary  mode  is  obtained  by 
substituting  equation  (8)  into  equation  (6)  and  integrating  to  obtain: 

Dyy  -  1  -  (k2c  ?fJ)  -  -  4up2  l  ANy/(u2  -  nV)  (1?) 

where: 

ANv  -  (N2/s2)exp(-s?/2)IN(s2/2)  +  (nt>/np)[(JN*  )2  +  (sxb/2)JNMNH] 

The  function  A^y  is  undulatory  about  zero,  therefore,  neetinq  the 
instability  condition  of  Tataronis  and  Crawford  (1970).  Using  a  four- 
term  approximation  to  the  dispersion  relation,  instabilities  are  found 


13 


for  (see  Appendix  R): 


nr>/np  =  TO  and  xb  =  9 
such  that: 

to/fi  -  1 .  *>  +0.2i 

and: 

nb/np  =  ^  and  xb  =  3 
such  that: 

i.j/n  =  1.8  +  n.fii 

Here  it  is  shown  that  unstable  solutions  exist  within  the  desired  fre¬ 
quency  ranqe.  The  larqer  beam-to-backqround  density  ratios  required 
to  support  the  extraordinary  mode  are  consistent  with  the  work  of 
Ohnuna  et  al .  (1981).  j 

C.  The  Ordinary  Mode  j 

The  dispersion  relation  for  the  ordinary  node  is  obtained  by  ! 

substituting  equation  (8)  into  equation  (7)  and  integrating  to  obtain: 

n2Z  =  1  -  (k2cW)  -  (Aoz  -  2wp2  l  ANz/(u»Z  -  hV)  (13) 

where: 

ANz  =  exp(-s2/2)IN(s2/2)  +  (nb/np)JN2 
Since  A^z  is  always  positive,  no  unstable  solutions  are  anticipated. 

CONCLUSIONS 

The  analytical  solutions  presented  here  show  that  electromagnetic 
non-resonant  instabilities  ma.y  be  excited  in  a  plasma  represented  as 
a  Maxwellian  background  with  a  monoenergetic  beam  of  electrons.  The 
electromagnetic  dispersion  relation  was  solved  for  kB  =  0,  using  the 
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approximation  <<  k*x  .  These  simplifications  permit  the  elements 
of  the  dispersion  relation  to  he  separated  and  solved  as  three  distinct 
eigennodes  of  perpendicularly  propaqatinq  waves: 

1.  Bernstein  Mode;  £  1  Ro,  k_  ii  F_ 

2.  Extraordinary  Mode;  £  l  Ro,  k  i  F 

3.  Ordinary  Mode;  £  n  Bo,  _k  j.  £ 

In  the  first  two  cases  it  is  shown  that  cyclotron  harmonic  modes 
nay  couple  between  the  first  two  harmonics  of  the  electron  gyrofreguen- 
cy,  excitinq  non-convecti ve  instabilities  with  qrowth  rates,  -y,  that 
are  a  finite  fraction  of  uir.  In  the  third  case  no  unstable  solutions 
exist.  Sample  solutions  presented  for  the  Rernstein  and  extraordinary 
nodes  are  shown  to  depend  on  the  ratios  n^/Pp  and  x^,  =  v^/vy.  These 
ratios  must  be  qreater  than  one. 

The  ratio,  nj,/np,  for  the  nearly  electrostatic  Bernstein  node 
can  be  compared  quantitatively  to  the  a  of  Tataronis  and  Crawford 
(1970)  and  the  ratio  Nc/N^  of  Ashour-Abdal 1  a  et  al.  (1975)  by  considerinn 
the  Maxwel 1 ian  portion  of  the  electron  distribution  as  the  "cold" 
component  and  the  delta-function  as  the  "hot"  or  "rinq"  component.  The 
results  for  n^/op  =  4.5  and  x^  =  10  agrees  with  the  results  for  a  =  0.2 
and  Nc/Nh  =  0.2.  These  are  values  associated  with  non-convective 
electrostatic  instabilities  at  frequencies  between  the  first  two  har¬ 
monics  of  the  electron  qyrof requency . 

For  the  extraordinary  mode  a  larger  density  ratio,  n^/np  «  10, 
is  required  to  excite  instabilities,  consistent  with  the  results  of 
Ohnuma  et  al.  (1981).  These  instabilities  may  be  excited  for  velocity 
ratios  of  3  and  9,  representing  ratios  of  the  "artificial"  and  "natural" 
electron  beam  densities,  respectively,  to  the  ambient  density.  These 
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electromagnetic  instabilities  are  non-convecti ve,  non-resonant,  with 
large  growth  rates.  They  may  be  excited  by  electron  beams,  given 
sufficient  bean-to-amhient  density  and  velocity  ratios.  The  ful  1 
range  of  parameters  over  which  they  may  occur  will  be  explored  through 
numerical  calculations.  For  further  understanding  of  the  rel ati onsln p 
of  these  waves  to  effects  observed  in  the  SCATHA  data  and  in  the  naqneto- 
-sphere,  the  full  electromagnetic  dispersion  relation,  including  the  kn 
terns,  must  be  solved  for  conditions  pertaining  to  space  vehicles  in 
space  plasmas. 
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APPENDIX  A 


Solution  for  Bernstein  Modes: 

The  dispersion  equation  for  the  Bernstein  modes  is: 

nxx  =  1  '  l  N^ANx/(u?  -  NV) 

where: 

aNx  =  s”2exp(-s2/2)IN(s2/2)  +  (nb/npsxb)JN(sxb)JN' (sxb) 
This  may  be  approximated: 

0XX  =  1  -  (4cob2/n2)[(A1«2/(u)2  -  SI2))  +  (4A?S12/(u,2  -  4S22))] 

Setting  Dxx  =  0  leads  t.o  the  fourth  order  equation: 

-  [5  +  (4cot)2/Sl2)(A1+4A?)]  u2/SJ2  + 

4[1  +  (4^2/rt2)  (Ax  +  A2)]  =  0 


( A  —  1 ) 


(A-2) 


(A— 3) 
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This  may  be  solved  as  a  quadratic,  then  transformed  to  polar  coordi nates, 
t  o  ol ita  i  n: 

m/a  =  (Kr)  l/2[cos(  0/2  +  nn)  +  isin(o/2  +  mi)],  n  =  0,1  (A-4) 

where: 

K  =  [S  +  (4<op2/.q?)(A1  +  4A2)]/2 
r2  =  x2  +  y2 
0  =  tan“l(.y/x) 
x  =  1 

V  =  [0(1  +  (4UpW)(Al  +  A2))/(S  +  (4^?AV>)(A1  +  4A?))2]  -  l]l/2 

For  instability,  the  followinq  condition  rust  be  satisfied: 

0  <  [5  +  (4^2/P2)(A1  +  4A?)]2  <  1 6  [1  +  (4Up2/P.2)  ( ^  +  A2)]  (A-S) 

A |  and  Ap  are  evaluated  as  functions  of  np/np  usinq  tabulated  values  of 
Bessel  functions  fron  Ahranowitz  and  Stequn  (1970)  for  the  approximate 
ranqe  of  the  arqument,  as  defined  by  the  constraints  on  sxp.  Substitution 
these  and  the  value  of  up/p  into  the  inequality  (A- 5),  determines  the 
constraints  on  np/np. 

For  sxp  =  3.6;  xp  =  20;  np/np  =  2: 

Ai  =  0.0464;  Ap  =  -0.0132. 

Substitutinq  in  equation  (A-4)  and  solving  for  o  and  Kr: 

•  1.95(0.89  +  0.45i)  =  1.73  +  0.886i  (A-6) 

A  second  solution  is  found  for  sxp  -  3.2,  xp  =  10.  Then  with 
np/nD  =  4.5;  A}  =  0.078,  Ap  =  -0.0223;  and 


id/ Q  *  1.83  +  1.14i 


APPENDIX  R 


Solution  for  Extraordinary  Modes: 

The  dispersion  equation  for  the  extraordinary  node  is: 

nyy  =  i  -  (kW)  -  \  ANy/(«?  -  nz2?-)  (R-i) 

where: 

ANy  =  (N?/s2)exp(-s?/2)lN,(s2/2)  +  (2nb/np)  [ ( J rJ ’  (sxb)  )2  + 
sxhJN’ (sxb)JM"(sxb)] 

This  nay  be  approximated: 

nyy  =  i  -  (k2cW)  -  - 

-  «2))  +  (A ?/(J  -  4ft2))]  (R-2) 

Settinq  nyy  =  0  leads  to  the  sixth  order  equation: 

U4  -  5 (wft)2  +  4n4)[(w2/k2c2)  -  1  -  [2^\/k?c2)l 

-  (4 c^j^/k 2C2) [«4( Aj+A2 )  -  M?(4A1+A2)]  =  n  (R-3) 

This  may  be  reduced  to  a  fourth  order  equation  by  usinq  the  approximation 
u?-/k2c?-  <<  1: 

[1  +  (2Wp2A0/k2c2)  +  (4<oD2/k2c2)(A1+A2)]  (u/ft)4  - 
[5(1  +  (2<^2A0/k2c2))  +  (4lop2/k2c2)(4A1+A2)]  (u>/ft)2  + 

4[1  +  (2  ui^hjk^c2)  ]  =  0  (R-4) 

This  nay  be  solved  as  a  quadratic,  then  transformed  to  polar  coordinates, 
analoqous  to  the  solution  for  the  Rernstein  mode  of  Appendix  A,  to  obtain 

uifsi  =  (Kr) l/2[cos( 0/2  +  mir)  ±  isin(o/2  +  mu)],  m  =  0,1  (R-5) 


where: 


-  |  ‘.(i  <  ^■^/A0A v,«/)  .  (A-^/kVj'AA j<a2).!/ 

?[1  +  ?c2)  *  (4ab?/k2c2)(A1+A2)] 

r2  =  x2  +  y2 
0  =  tan-l(y/x) 
x  =  l 

y  -  [[16(1  +  2Wp2A0/k2c2)[l  +  (?Wp?A0/k2c2)  +  (4Wp2/k2c2)(A1+A2)]/ 
[6(1  +  ?,^2A0/k2c2)  +  (4,^2/k2c2)(4A1+A?)]2]  -  I]1/2 

For  instability,  the  followinq  condition  that  must  be  satisfied: 

0  <  [6(1  +  2wp2AQ/k2c2)  +  (4<up2/k2c2)(4A1+A;?)]2 

<  16(1  ♦  ?(Jp2A0/k2c2)[l  +  ?(l)p2A0/k2c2  +  (4n,2/k2c2)(A1+A?)]  (6-6) 

An  and  An+|  are  evaluated  as  functions  of  n^/np  usinq  the  tabulated  val¬ 
ues  of  Vessel  functions  from  Abranowitz  and  Steaun  (1970)  for  the 
appropriate  ranqe  of  the  argument  as  defined  by  the  constraints  on 
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sxb.  Substitute  these  and  the  value  of  uu  /k'c  into  the 
inequality  (6-6)  to  obtain  the  constraints  on  nfo/np. 

For  sxb  =  4.0;  x^  -  9.0;  n^/np  =  10: 

A0  -■  1.06,  A|  =  -0.714,  A2  =  2.72.  Substitutinq  and  solvinq  for 
0  and  Kr; 

u/.0  =  1.66(0.99  +  0. 1 4  i )  =  1.63  +  0.22i  (R-7) 

A  second  solution  is  found  for  sx^  =  4.4;  xi-,  =  3.  Then  with 
nb/np  =  40  =  3.06,  A^  =  -2.4,  Ap  =  2.79,  and 


u/H  -  1.8  +  D.6i 


